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Abstract 

Every  economic  model  should  include  an  estimate  of  its  stability  and  predictability.  A  new  measure,  the  first  passage 
time  (FPT)  which  is  defined  as  the  time  period  when  the  model  error  first  exceeds  a  pre-determined  criterion  (i.e.,  the 
tolerance  level),  is  proposed  here  to  estimate  the  model  predictability.  A  theoretical  framework  is  developed  to  deter¬ 
mine  the  mean  and  variance  of  FPT.  The  classical  Kaldor  model  is  taken  as  an  example  to  show  the  robustness  of  using 
FPT  as  a  quantitative  measure  for  identifying  the  model  stability. 
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1.  Introduction 

One  of  the  most  interesting  theories  of  business  cycles  in  the  Keynesian  vein  is  that  expounded  in  a  pioneering  article 
by  Kaldor  [1].  It  is  distinguishable  from  most  other  contemporary  treatments  since  it  utilizes  non-linear  functions, 
which  produce  endogenous  cycles,  rather  than  the  linear  multiplier-accelerator  kind  which  rely  largely  on  exogenous 
factors  to  maintain  regular  cycles. 

The  savings  and  investment  functions  in  the  income-expenditure  theory  of  Keynes  are  linear.  Kaldor  [1]  proposed 
that  the  treatment  of  savings  and  investment  as  linear  curves  simply  does  not  correspond  to  empirical  reality.  Let  ( Y,K ) 
denote  gross  product  and  capital  stock.  The  investment  function  I(Y,K)  and  savings  function  S(Y,K)  are  increasing 
functions  with  respect  to  Y.  Over  the  savings  and  investment  functions,  Kaldor  superimposed  Keynes’s  multiplier  the¬ 
ory,  namely,  that  gross  product  changes  to  clear  the  goods  market. 

Chang  and  Smyth  [2]  and  Varian  [3]  translated  Kaldor’s  trade  cycle  model  into  more  rigorous  context:  the  former 
into  a  limit  cycle  and  the  latter  into  catastrophe  theory.  Output,  as  we  saw  via  the  theory  of  the  multiplier,  responds  to 
the  difference  between  savings  and  investment.  If  there  is  excess  goods  demand  (which  translates  to  saying  that  invest¬ 
ment  exceeds  savings,  I  >  S),  then  gross  product  rises  ( dYIdt  >  0),  whereas  if  there  is  excess  goods  supply  (which  trans¬ 
lates  to  savings  exceeding  investment,  I  <  S),  then  Y  falls.  The  Kaldor  system  is  represented  by 
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where  n  is  the  adjustment  coefficient  in  the  goods  market. 

It  is  widely  recognized  that  the  uncertainty  can  be  traced  back  to  three  factors:  (a)  measurement  errors,  (b)  model 
errors  such  as  uncertain  model  parameters,  and  (c)  chaotic  dynamics.  Measurement  errors  cause  uncertainty  in  initial 
conditions.  Discretization  causes  truncation  errors.  The  chaotic  dynamics  may  trigger  a  subsequent  amplification  of 
small  errors  through  a  complex  response. 

Traditionally,  the  small  amplitude  stability  analysis  (linear  error  dynamics)  is  used  to  study  the  model  stability.  This 
method  is  divided  into  two  steps.  First  step  is  to  find  equilibrium  states  of  the  dynamic  system.  The  second  step  is  to 
investigate  temporal  evolution  of  perturbations  relative  to  the  equilibrium  states.  It  is  well  known  that  comparative  sta¬ 
tic  analysis  is  only  valid  if  there  is  a  tendency  of  the  variables  towards  the  new  equilibrium.  If  such  movement  does  not 
take  place,  the  comparative  static  exercises  do  not  really  contribute  to  the  knowledge  of  the  evolution  of  the  economic 
systems  [9],  In  modeling  the  endogenous  business  cycles,  the  macroeconomic  equilibrium  is  never  reached,  but  oscilla¬ 
tory  motion  is  usually  found  around  it.  Therefore,  it  is  more  appropriate  to  examine  the  comparative  dynamics  of  the 
systems  (la)  and  (lb),  i.e.,  to  investigate  the  whole  change  over  time  caused  by  a  change  in  initial  conditions,  exogenous 
parameters  or  reaction  coefficients  [4],  In  the  context  of  business  cycle  theory,  the  change  of  amplitude  of  business  cycle 
is  caused  by  the  changes  in  the  parameters  and  initial  conditions. 

A  question  arises:  how  long  is  the  Kaldor  model  (1)  valid  since  being  integrated  from  its  initial  state?  This  has  great 
practical  significance.  For  example,  if  the  model  validity  time  shorter  than  the  business  cycle,  the  model  does  not  have 
any  capability  to  predict  the  business  cycle.  In  this  paper,  probabilistic  stability  analysis  is  proposed  to  investigate  the 
model  valid  period.  This  method  is  on  the  base  of  the  first-passage  time  (FPT)  for  model  prediction. 


2.  FPT  for  prediction 

Let  an  iV-dimensional  vector,  x(t)  =  [x(l)(t),  x(2Hf), • .  .,x(iV)(0]  represent  a  set  of  economic  variables  governed  by 
dx 

^  =  (2) 

which  is  the  extension  of  (la)  and  (lb).  Flere  f  is  a  functional.  Individual  economic  prediction  is  to  find  the  solution  of 
(2)  with  an  initial  condition 


x(t0)  =  x0.  (3) 

Uncertainty  in  economic  models  leads  to  the  addition  of  stochastic  forcing.  For  simplicity,  a  stochastic  forcing  (f')  is 
assumed  to  be  white  multiplicative  or  additive  noise,  and  (2)  becomes 

dx 

—  =  f(x,r) +f'(x,t),  f(x,t)  =  k(x,t)g(t),  (4) 

where  k(x,t)  and  g (t)  are  the  forcing  covariance  matrix  {k,-,-}  (dimension  of  Nx  N)  and  the  vector  delta-correlated  pro¬ 
cess  (dimension  of  N),  respectively. 

Let  x(f)  be  the  reference  solution  which  satisfies  (3)  with  the  initial  condition,  x(t0)  =  x0.  The  model  error  z  is  deter¬ 
mined  as 

Z  (t)  =x(t)-x'(t),  (5) 

where  x'(?)  is  one  of  individual  prediction.  Two  vectors  x(f)  and  x’(t)  are  considered  as  reference  and  prediction  points 
in  the  iV-dimensional  phase  space. 

To  quantify  FPT,  we  first  define  two  model  error  limits.  First,  the  forecast  error  cannot  be  less  than  a  minimum  scale 
d,  which  depends  on  the  intrinsic  noises  existing  in  the  model.  Second,  the  forecast  error  cannot  be  more  than  a  max¬ 
imum  scale  (tolerance  level)s.  The  prediction  is  valid  if  the  reference  point  x(f)  is  situated  inside  the  ellipsoid  ( Se ,  called 
tolerance  ellipsoid)  with  center  at  x'(t)  and  size  £.  When  x(f)  coincides  with  x'(/),  the  model  has  perfect  prediction.  The 
prediction  is  invalid  if  the  reference  point  x(f)  touches  the  boundary  of  the  tolerance  ellipsoid  at  the  first  time  from  the 
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Fig.  1 .  Phase  space  trajectories  of  model  prediction  y  (solid  curve)  and  reality  x  (dashed  curve)  and  error  ellipsoid  SE(t)  centered  at  y. 
The  positions  of  reality  and  prediction  trajectories  at  time  instance  are  denoted  by  and  “O”  respectively.  A  valid  prediction  is 

represented  by  a  time  period  (t  —  t0)  at  which  the  error  first  goes  out  of  the  ellipsoid  Se(t). 


initial  state  that  is  FPT  for  prediction  (Fig.  1 ).  FPT  is  a  random  variable  when  the  model  has  stochastic  forcing  or  initial 
condition  has  random  error.  Its  statistics  such  as  the  probability  density  function,  mean  and  variance  can  represent  how 
long  the  model  can  predict.  The  FPT,  x  —  t  —  t0,  depends  on  the  initial  model  error,  z0  =  z(/0),  tolerance  level  e,  and 
model  parameters.  The  longer  the  FPT,  the  more  stable  of  the  economic  model  is. 


3.  Backward  Fokker-Planck  equation 

The  conditional  probability  density  function  (PDF)  of  FPT  with  a  given  initial  error, P[(t  —  t0)lzo]>  satisfies  the  back¬ 
ward  Fokker-Planck  equation  [5,6] 


dP  _  y-  /  ,  _  dxA  dP  _  1 
dt  fY  V  dt  J  Sz°  2 


N  N  N  ~2p 

/=i  /=  l  /=  l  1  j 


=  0, 


(6) 


where  the  coefficients  kg  are  the  components  of  the  forcing  covariance  matrix  k(x,  t)  and  (z^z®, . .  ,z°N)  are  the  compo¬ 
nents  of  the  initial  error  z0.  Integration  of  PDF  over  t  leads  to 


P[(t  -  fo)|zo]df 


=  1. 


(7) 


The  A'th  FPT  moment  {k-  1,2,  . . .)  is  calculated  by 


P  oo 

xk(z0 )  =  k  P[(t-t0)lzo](f-fo)t"1df,  k  =  1, . . . ,  oo. 

JtQ 


(8) 


If  the  initial  error  z0  reaches  the  tolerance  level,  the  model  loses  prediction  capability  initially  (i.e.,  FPT  is  zero) 
P[(t- ?o)|zo]  =  0  at  J(z0)  =  e2,  (9a) 
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which  is  the  absorbing  type  boundary  condition.  Here  /( z0)  denotes  the  norm  of  z0.  If  the  initial  error  reaches  the  noise 
level  the  boundary  condition  becomes  [6] 

0P[(t-  f0)|zo] 


0c 


!J) 


=  0  at  J(zo)  =  <5  , 


(9b) 


which  is  the  reflecting  boundary  conditions.  Here,  £  is  the  noise  level.  Usually, 

S  <  E. 

Mean,  variance,  skewness,  and  kurtosis  of  the  FPT  are  calculated  front  the  first  four  moments 


(t)  = 

(10a) 

{St2)  =  t2  -  t2, 

(10b) 

(<5t3)  =  r3  —  3t2ti  +  2tj, 

(10c) 

(0T4)  =  r4  —  4t3T!  +  6t2t2  —  3t[, 

(10d) 

where  the  bracket  denotes  the  ensemble  average  over  realizations  generated  by  stochastic  forcing. 


4.  Kaldor  model 


4.1.  Average  method 


The  Kaldor  model  ( 1 )  can  be  written  by 


d 

Y 

=  A 

Y' 

+ 

MY,Ky 

d  t 

K 

K 

MY,K)_ 

(11) 


where 


A  = 


'au 

a\2 

_a2 1 

a22  _ 

H(IY  —  SY)  p{Ik  —  SK ) 
h  Ik 


J  (V,K') 


(12) 


is  the  Jacobian  matrix  evaluated  at  the  equilibrium  (Y*,K*)  =  (0,0),  and  f\  and  f2  are  nonlinear  terms.  Here,  the  sub¬ 
scripts  denote  the  partial  differentiation;  —IK  is  the  depreciation  coefficient;  IY  is  the  coefficient  for  reinvested  profits. 
The  Jacobian  matrix  A  has  a  determinant: 


|A|  —  p(IY  —  SY)1K  —  p(IK  —  SK)IY  —  n(SKIr  —  IkSy),  (13) 

where,  since  1K  <  0  and  SK,  SY,  /y>  0  then  |A|  >  0,  thus  we  have  regular  (non-saddle  point)  dynamics. 

Following  Chiarella’s  [7]  approach  of  using  the  polar  coordinates,  Y  =  r  cos  6,  K  =  rsin0,  the  Kaldor  model  is 
transformed  into 


d;- 
d  t 
dd 
d  t 


=  ;-[anCOS20  +  a22sin  9  +  (a12  +  a2\)  cos  0  sin  0]  +  gx(r,  8)  cos  9  +  g2(r,  8)  sin  9, 

-2  n  ,  \  ■  n  n  \s2  (ri  9)  cos  9  —  g,  (r,  9)  sin  01 

=  a2icos20  —  fli2sm“0  +  (a22  —  an)  sm0cos  0  +  - 1 - - - - , 


where 


(14a) 

(14b) 


gi(r,9)  =  /)(rcos0,rsin0),  i  =  1,2. 

The  right-hand  sides  of  Eqs.  ( 14a, b)  are  periodic  in  0  and  can  be  expanded  into  Fourier  series  with  desired  degree  of 
accuracy.  The  first  term  of  this  expansion  is  obtained  by  averaging  the  right-hand  sides  of  Eqs.  (14a)  and  (14b)  with 
respect  to  0  on  the  interval  [0,2rc].  In  this  way,  the  equation  for  averaged  amplitude  of  the  business  cycle  [7] 


dp  _  Tr(A) 

d7=  2 


p  +  G(p),  Tr(A)  —  p(Ir  —  SY)  +  IK  >  0, 


(15) 
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where  p  is  the  first-order  approximation  to  r  and 

G(p)=2jt  J  [S\(P,S)cose  +  g2(p,0)smd}de.  (16) 

The  nonlinear  term  G (p)  may  have  several  different  forms.  Here,  we  use  the  profit-capital-accumulation  dynamic  model 
as  an  example  for  illustration.  Semmler  [8]  developed  an  endogenous  cycle  model  of  profit  and  capital  accumulation 
dynamics  on  the  base  of  the  Kaldor  dynamics.  The  Kaldor  nonlinear  term  G(p)  has  the  following  form  [9]: 

G(p)  =  -^p\  (17) 

The  coefficients  (p,  IY ,  Ik,  Sy)  in  Eq.  (15)  are  usually  difficult  to  measure.  This  may  lead  to  the  addition  of  stochastic 
forcing  to  the  business  cycle  model.  For  simplicity,  the  stochastic  forcing  is  assumed  to  be  white  multiplicative  or  addi¬ 
tive  noise,  and  Eq.  ( 1 5)  becomes 

^  =  op  +  G(p)  +  m(t)p,  a  =  Tr^  ,  (18a) 

where  m(t)  is  a  random  parameter  with  zero  mean  and  pulse-type  variance, 

( m(t ))  =  0,  (. m(t)m(t '))  =  q2A(t  —  t1),  (18b) 

where  the  bracket  indicates  the  ensemble  average;  q  represents  the  strength  of  the  uncertainty  in  coefficients  in  a  (such  as 
in  the  depreciation  coefficient  and  the  coefficient  for  reinvested  profits);  and  A (t)  is  a  delta  function. 


4.2.  FPT  moments 

Let  p(t)  be  the  reference  solution  which  satisfies  (18a)  with  the  initial  condition,  p(t0)  -  p0.  The  forecast  error  Z  is 
determined  as 

Z(t)  ~p{t)-p'{t), 

where  p'(t)  is  one  of  individual  prediction  corresponded  to  perturbing  initial  condition  and/or  stochastic  forcing.  The 
backward  Fokker-Planck  Eq.  (6)  for  this  case  is  simplified  to 


0P 

dt 


dP 


1  ^ 2p  n 

-q  — T  =  0. 

>  1  rsyl 
^0 


(19) 


We  multiply  Eq.  (19)  by  (t  —  /0)n,  integrate  with  respect  to  1  from  t0  to  oo,  use  the  condition  (7),  and  obtain  the  equa¬ 
tions  of  the  «th  moment  of  FPT 


2  dZl 


+  [aZo  -fzj] 


P  73) 

dZ0 


=  T0  —  1 . 


(20) 


Eq.  (20)  is  linear,  time-independent,  and  second-order  differential  equations  with  the  initial  error  Z0  as  the  only  inde¬ 
pendent  variable.  Two  boundary  conditions  for  ti  and  r2  can  be  derived  from  (9a)  and  (9b) 


rk  =  0,  for  |Z0|  =  s; 

=  0.  for  |Z0|  =  S. 
0ZO  1  1 


Analytical  solution  of  (20)  with  the  boundary  conditions  (21a)  and  (21b)  is 

t„{zo,i,E,a,q2)  =4  J  y exp  j  J r  nz„-i{x)x?~- Qxp^-^x^jdx 


d y, 


(21a) 

(21b) 

(22) 
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where  i0  =  1,  and 


(23) 


are  non-dimensional  initial  error  and  noise  level  scaled  by  the  tolerance  level  e.  The  moments  of  FPT  depend  on  two 
types  of  parameters:  (a)  prediction  parameters  (z0,<^,e),  and  (b)  model  parameters  (a-,q2,fi),  It  is  convenient  to  use  the 

two  lowest  order  statistics,  mean  FTP  ( )  and  the  standard  deviation  of  FTP  |s  =  sj {<5t2)J ,  to  analyze  the  stability  of 
the  Kaldor  model. 


5.  Model  predictability 

The  first  two  moments  z1  and  r2  can  be  taken  as  the  stability  measure  of  the  dynamic  system.  The  longer  the  mean 
FPT,  the  more  the  stable  of  the  system  is.  The  dependence  of  tyand  r2  on  the  two  types  of  parameters  is  investigated 
separately.  The  three  model  parameters  are  taken  as  cf  =  0.2,  a  =  1.0,  n  =  1. 


Fig.  2.  Contour  plots  of  versus  (r0,  {)  for  four  different  values  of  e  (0.1. 0.2. 1.0,  2.0)  using  the  Kaldor  model  with  given  model 

parameters  a  =  1.0,  q  =  0.2,  pi  =  1.  The  contour  plot  covers  the  half  domain  due  to  ro  > 


P  C.  Chu  /  Chaos,  Solitons  and  Fractals  27  (2006)  1355-1368 


1361 


Fig.  3.  Contour  plots  of  t2(zo>  {,  £)  versus  (z0,  £)  for  four  different  values  of  e  (0. 1 , 0.2, 1 .0, 2.0)  using  the  Kaldor  model  with  given  model 
parameters  a  =  1.0,  q2  =  0.2,  /( =  1.  The  contour  plot  covers  the  half  domain  due  to  zq  >  {■ 


Figs.  2  and  3  show  the  contour  plots  of  T\(z0,  £,e)  and  x2(z0,£,e)  versus  (z0,<£)  for  four  different  values  of  e  (0.1,  0.2, 
1.0,  and  2.0).  Following  features  can  be  obtained:  (a)  for  given  values  of  (z0,£)  [i-e.,  the  same  location  in  the  contour 
plots],  both  Tj  and  r2  increase  with  the  tolerance-level  e.  (b)  For  a  given  value  of  tolerance-level  e,  both  x\  and  r2  are 
almost  independent  on  the  noise  level  £  (contours  are  almost  paralleling  to  the  horizontal  axis)  when  the  initial  error  (z0) 
is  much  larger  than  the  noise  level  (£).  This  indicates  that  the  effect  of  the  noise  level  (<f)  on  xt  and  i2  becomes  evident 
only  when  the  initial  error  (z0)  is  close  to  the  noise  level  (£).  (c)  For  given  values  of  (e,  £),  both  x\  and  x2  decrease  with 
increasing  initial  error  z0. 

Figs.  4  and  5  show  the  curve  plots  of  t^zo,  £,e)  and  s(z0,£,e)  versus  z0  for  four  different  values  of  tolerance  level,  s 
(0.1,  1,  2,  and  3)  and  four  different  values  of  random  noise  £  (0.1,  0.2,  0.4,  and  0.6).  Following  features  are  obtained:  (a) 
x i  and  .v  decrease  with  increasing  z0,  which  implies  that  the  higher  the  initial  error,  the  shorter  the  mean  FPT  (or  lower 
model  predictability)  and  the  smaller  the  s  (or  lower  variability  of  the  model  predictability)  are;  (b)  x\  and  s  decrease 
with  increasing  noise  levels,  which  implies  that  the  higher  the  noise  level,  the  lower  the  FPT  and  its  variability  are;  and 
(c)  t i  and  s  increase  with  increasing  e,  which  implies  that  the  higher  the  tolerance  level,  the  longer  the  FPT  (or  higher 


1362 


PC.  Chu  /  Chaos.  Solitons  and  Fractals  27  (2006)  1355-1368 


Fig.  4.  Dependence  of  t, (z0,  (,  s)  on  the  initial  condition  error  z0  for  four  different  values  of  e  (0.1, 1.0, 2.0, 3.0)  and  four  different  values 
of  random  noise  £  (0.1, 0.2, 0.4, 0.6)  using  the  Kaldor  model  with  given  model  parameters  a  =  1.0,  q2  =  0.2,  /( =  1. 


model  predictability)  and  the  larger  the  s  (or  higher  variability  of  the  model  predictability).  For  the  same  values  of 
model  parameters,  the  standard  deviation  of  FTP  is  much  smaller  than  the  mean  FTP.  Note  that  the  results  presented 
in  this  subsection  is  for  a  given  values  of  growth  rate  (a  =  1.0),  stochastic  forcing  ( q 2  =  0.2),  and  the  Kaldor  model 
parameter  (p  =  1)  only. 


6.  Dependence  of  t|  and  s  on  the  Kaldor  model  parameters 

6.1.  Dependence  on  the  growth  rate  o 

To  investigate  the  sensitivity  of  %\  and  s  to  the  model  parameter  a,  the  other  two  model  parameters  are  kept 
unchanged  (q2  —  0.2,  p=  1).  The  model  parameter  a  takes  values  of  0.25,  0.5,  and  1.0.  Figs.  6  and  7  show  the 
curve  plots  of  Ti (z0,^,o)  and  s(z0,i,ct)  versus  z0  for  two  tolerance  levels  (e  =  0.05,0.25),  two  noise  levels 
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Fig.  5.  Dependence  of  ,s(r0,  e)  on  the  initial  condition  error  z0  for  four  different  values  of  £  (0. 1 , 1 .0, 2.0, 3.0)  and  four  different  values 
of  random  noise  £  (0.1, 0.2, 0.4, 0.6)  using  the  Kaldor  model  with  model  parameters  a  =  1.0,  q2  =  0.2,  g  =  1. 


(£  =  0.1, 0.6),  and  three  different  values  of  a  (0.25,  0.5,  and  1.0).  tj  and  s  decrease  with  increasing  a  for  all  com¬ 
binations  of  noise  level  (£  =  0.1, 0.6)  and  tolerance  level  (s  —  0.05,0.25).  This  indicates  that  increase  of  Tr(A)  [i.e., 
tr]  destabilizes  the  Kaldor  model  (decreasing  r()  and  decreases  the  variability  of  the  model  predictability  (decreas¬ 
ing  s);  and  decrease  of  Tr(A)  stabilizes  the  Kaldor  model  (increasing  t])  and  increases  the  variability  of  the  model 
predictability  (increasing  5). 


6.2.  Dependence  on  the  stochastic  forcing  cf 

To  investigate  the  sensitivity  of  and  s  to  stochastic  forcing  q1  and  the  Kaldor  model  parameter  ft,  the  growth  rate 
is  kept  unchanged  (a  =  1).  Figs.  8  and  9  show  the  curve  plots  of  (z0,  £,  q  )  and  s(z0,^,q2)  versus  z0  for  two  adjustment 
coefficients  (g  =  0.01,  10),  two  noise  levels  (£  =  0.1,  0.6),  and  three  different  values  of  q  (0.1,  0.25,  and  0.5)  representing 
weak,  normal,  and  strong  stochastic  forcing.  Two  regimes  are  found:  (a)  zi  and  s  decrease  with  increasing  cf  for  large 
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Fig.  6.  Dependence  of  t,(z0 ,^,q2)  on  the  initial  condition  error  r0  for  three  different  values  of  a  (0.25,0.5, 1.0)  and  given  values  of  cf 
(=0.2)  and  /t(=l)  using  the  Kaldor  model  with  two  different  values  of  e  (0.05,0.25)  and  two  different  values  of  noise  level  {  (0.1, 0.6). 
Note  that  the  mean  FTP  is  not  sensitive  to  a. 


noise  level  (£  =  0.6),  (b)  and  s  increase  with  increasing  cf  for  small  noise  level  (£  =  0.1)  and  (c)  both  relationships 
(increase  and  decrease  of  z\  and  s  with  increasing  q2)  are  independent  of  /(. 


7.  Stabilizing  and  destabilizing  regimes 


This  indicates  the  existence  of  stabilizing  and  destabilizing  regimes  of  the  dynamical  system  depending  on  stochastic 
forcing.  For  a  small  noise  level,  the  stochastic  forcing  stabilizes  the  dynamical  system  and  increase  the  mean  FPT.  For  a 
large  noise  level,  the  stochastic  forcing  destabilizes  the  dynamical  system  and  decreases  the  mean  FPT.  The  two  regimes 
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Fig.  7.  Dependence  of  ,s(r0,  i,q2)  on  the  initial  condition  error  r0  for  three  different  values  of  a  (0.25,0.5, 1.0)  and  given  values  of  q1 
(=0.2)  and  ^(=1)  using  the  Kaldor  model  with  two  different  values  of  e  (0.05,0.25)  and  two  different  values  of  noise  level  {  (0.1, 0.6). 
Note  that  the  standard  deviation  of  FTP  is  not  sensitive  to  a. 


can  be  identified  analytically  for  small  tolerance  level  ( e  — »  0).  The  initial  error  z0  should  also  be  small  (z0  ~  e).  The  solu¬ 
tions  (22)  becomes 


LimT^zo,^,  e,  tr,^2)  = 


1 


In 


a  -  q2/2  I  \z0J  2a  —  q2 


T  r#”‘ 


-  1 


(24) 


The  Lyapunov  exponent  is  identified  as  ( a  —  q2 / 2)  for  dynamical  system  (16)  [4],  For  a  small  noise  level  (c,  1),  the 
second  term  in  the  bracket  of  the  right-hand  of  (24) 


2a  —  q 


jf-t 


-  1 


R  -  - 


(25) 
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Fig.  8.  Dependence  of  ii(r0 .f,?2)  on  the  initial  condition  error  z0  for  three  different  values  of  q 1  (0.1,0.25,0.5)  and  given  parameters 
(a  =  1 .0,  £  —  0.2)  using  the  Kaldor  model  with  two  different  values  of  the  adjustment  coefficient  (0.01, 10)  and  two  different  values  of 
noise  level  £  (0.1, 0.6).  Note  that  the  mean  FTP  is  not  sensitive  to  varying  fi. 


is  negligible.  The  solution  (24)  becomes 


Lim  Ti(z0,  £,£,  a,q2) 
£—►0 


(26) 


which  shows  that  the  stochastic  forcing  (q  #  0),  reduces  the  Lyapunov  exponent  (a  —  q2! 2),  stabilizes  the 
dynamical  system  (15),  and  in  turn  increases  the  mean  FPT.  On  the  other  hand,  the  initial  error  z0  reduces  the  mean 
FPT. 

For  a  large  noise  level  f,  the  second  term  in  the  bracket  of  the  right-hand  of  (24)  is  not  negligible.  For  a  positive 
Lyapunov  exponent,  2er  —  q1  >  0,  this  term  is  always  negative  [see  (25)].  The  absolute  value  of  R  increases  with  increas¬ 
ing  q2  (remember  that  £  <  c0  <  1).  Thus,  the  term  ( R )  destabilizes  the  stochastic  Kaldor  system  (15),  and  reduces  the 
mean  FPT. 
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Fig.  9.  Dependence  of  s(z0,  q1)  on  the  initial  condition  error  r0  for  three  different  values  of  q 2  (0.1,0.25,0.5)  and  given  parameters 

(a  =  1.0,  e  =  0.2)  using  the  Kaldor  model  with  two  different  values  of  the  adjustment  coefficient  /t  (0.01, 10)  and  two  different  values  of 
noise  level  £  (0.1, 0.6).  Note  that  the  mean  FTP  is  not  sensitive  to  varying  /t. 


8.  Conclusions 

(1)  The  Kaldor  model  stability  and  predictability  are  not  only  affected  by  the  model  parameters  such  as  the  depre¬ 
ciation  coefficient,  the  coefficient  for  reinvested  profits  [i.e.,  Tr(A)],  the  adjustment,  but  also  by  the  prediction 
parameters  (such  as  initial  error,  tolerance  level,  and  noise  level.  The  capability  of  the  FPT  approach  in  evalu¬ 
ating  model  stability  and  predictability  is  demonstrated  using  the  nonlinear  Kaldor  model. 

(2)  Uncertainty  in  economic  models  is  caused  by  uncertain  measurements,  computational  accuracy,  and  uncertain 
model  parameters.  This  motivates  to  the  inclusion  of  stochastic  forcing  in  economic  models  such  as  the  Kaldor 
model.  The  backward  Fokker-Planck  equation  can  be  used  for  evaluation  of  economic  model  stability  and  pre¬ 
dictability  through  the  FPT  calculation. 

(3)  A  theoretical  framework  was  developed  in  this  study  to  determine  various  FPT  moments  (t*),  which  satisfy  time- 
independent  second-order  linear  differential  equations  with  given  boundary  conditions.  This  is  a  well-posed  prob¬ 
lem  and  the  solutions  are  easily  obtained. 
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(4)  For  the  Kaldor  model,  the  following  features  are  detected  from  the  FPT  calculation:  (a)  decrease  of  ii  and  s  with 
increasing  initial  condition  error  (or  with  increasing  random  noise),  (b)  slow  increase  of  tx  and  s  with  increasing 
tolerance  level  e,  (c)  For  the  same  values  of  model  parameters,  the  standard  deviation  of  FTP  is  much  smaller 
than  the  mean  FTP,  which  indicates  that  the  first  moment  is  a  reliable  indicator  of  the  model  stability. 

(5)  Both  stabilizing  and  destabilizing  regimes  are  found  in  the  Kaldor  model  depending  on  stochastic  forcing.  For  a 
small  noise  level,  the  stochastic  forcing  stabilizes  the  Kaldor  system  and  increases  the  mean  FPT.  For  a  large 
noise  level,  the  stochastic  forcing  destabilizes  the  Kaldor  system  and  decreases  the  mean  FPT. 

(6)  Model  stability  depends  on  Tr(A).  Increase  of  Tr(A)  destabilizes  the  Kaldor  model  (decreasing  ti)  and  decrease  of 
Tr(A)  stabilizes  the  Kaldor  model  (increasing  ti).  The  model  stability  does  not  depend  on  the  Kaldor  model 
parameter  /r. 
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